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ABSTRACT 


The dangers of copyright protection can impact 2D vector maps, having a 
knock-on effect on the use of vector data. To achieve invariance property, 
uniform RST (rotation, scaling and translation) and disguising the digital 
vector map’s information by implementing distortion control, is all done by 
using watermarking schemes. Convert an original map, then engrain 
the watermark. An LCA algorithm is used in this study, as a newly proposed 
way to protect the vector maps under copyright. The procedure is operated in 
this order: 1) use an original map, altered by the LCA algorithm, 
2) use the coefficient of the transformation to engrain the watermark, 
inserting the resulting frequency into the LSB wave, 3) the watermarked map 
is acquired by using the inverse LCA map transformation. Further 
investigations discovered that the necessary standards of fidelity and 
invisibility can be achieved using this technique. This procedure also gives 
out numerous frequency domains for digital watermarking; as well as being 
resilient to signal and geometric invasions. 
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1. INTRODUCTION 

From the making of paper maps to digital data types, we have seen the progression of geospatial 
data over the past decades. Improvements in computer hardware relative to geographic data collection tools 
have caused this. Its proficiencies include GPS (geographic positioning systems) with satellites able to 
retrieve trusted spatial positioning data. The functions of analog data or print has been replaced by vector 
maps as the standard data units in GIS (geographic information systems) [1]. An essential component of GIS 
are vector maps, a type of GIS spatial data with many distinctive qualities, covers many areas from 
navigation to traffic data. Current investigations intend to copyright protect the owners of these systems, as 
they become more popular and used across the board. In comparison to multimedia data which has fixed 
relative positions, graphics have several independent components. Graphics can also make copying scarce 
due to its abundance of topology and engineering information. The incentive for research into DW (digital 
watermarking) for digital vector maps, can be formed from combining these factors [2]. 

Deciding on who owns the digital map and whether it has validity, is determined by the digital map 
stakeholders, due to DW. Alteration when the media point is altered and protection against data extraction are 
the robust benefits of DW. Robust DW for map vectors can be used on both the transformation and spatial 
domains, and a strong DW application can be classed as copyright protection [3]. Discrete wavelet transform 
(DWT) [4], discrete cosine transform (DCT) [5], and fast fourier transform (FFT) [6] are all included in 
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the principal transform algorithms. The convenient application of the spatial domain is at 
the expense of robustness and invisibility, therefore to copyright applications, robust DW procedures 
concentrate on the transformation domain [7-9]. 

Research on robust DW algorithms in the transformation domain is focused on images [10-15] and 
audio [16-20]. Vector maps differ from multimedia data because they are not limited to functions with 
integers, however different methods are needed to enter data into a vector map. IWT (integer wave 
transform) [21], DFT [2] and DCT [22, 23] can be used by robust DW in the vector map transformation (as 
proof has shown), therefore our new procedure is compared to the DWT, DCT and FFT techniques in 
the current investigation. For each co-ordinate sequence, the final 12 digits of the decimal are taken and then 
put into a matrix: (a) then altered into 4 bands, (b) different vector data methods associated with copyright 
digital maps have been examined, [24-26] including changeable watermarking in the spatial domain and [27] 
watermarking in the spatial topology domain. The progression of copyright protecting vector maps is the aim 
of this research. Implementation of the LCA transformation algorithm is the advanced technique for 
watermarking vector maps, as suggested by this paper. Multimedia watermarking has been using cellular 
automata transform (CAT) for a very long time [28-30]. Embedded media still has no research results in 
connection to its use on vector maps. Using the LCA algorithm has the benefits of high fidelity, great 
insertion results and high invisibility [31]. Data verification traits are present due to the use of the scrambling 
method [25], plus it is not limited to one transformation plane, making it unique. This sets it above 
the current frequency watermarking methods and it is able to accelerate DW via its multi-frequency domains. 

Our methods for the performance measurements used were: 1) NC calculation was used to measure 
the performance technique, 2) evaluation of quality, based on invisibility with RMSE calculations, 
3) fidelity with the longest distance, 4) Resistance and NC calculation against geometric attacks (e.g. RST). 
The new method functioned well, as is evident in the results of the invisibility test, which produced map test 
data RMSE values. NC data and distance change were within the accepted threshold, so map fidelity was 
high. There was a high level of resilience against geometric attacks, especially for the watermark as integrity 
was essential. There are four extra sections at the end of this paper. An explanation of the investigation 
procedures used in this paper are documented in section 2, followed by the results in section 3, then section 4 
has a conclusion of the study. 


2. RESEARCH METHOD 
2.1. Linear cellular automata 

A dynamical in a discontinuous frequency and time frame, symbolized by the LCA transform is 
shown below in (1). Each cell has a restricted group of states and form a lattice structure. Effectual and fast 
calculations of the discrete transformation are enabled by the algorithm [31]. 


(C'+?)T = Ma. (C‘)"(mod 2) (1) 


M,, represents the local transition matrix, given that n = 5 k: 


EELO O 0 0 0 
1111 0.. 0 0 0 
i ent he Lames Beene Bere 0 0 0 
0 1 11 1.. 0 0 0 
M, = 0 0 11 1.. 0 0 0 
O00 0 E O E 11 
00000.1111 
00000.1111 


if the transition matrix of a cellular automation (An) was represented by M,, the non-zero coefficients will 
become 1 as An is now the nth order penta-diagonal matrix. (Ct)? on the other hand, is a transposition of a 
linear matrix, composed of inter-changing and random binary numbers are shown in (2). 


(Cr) " = Mī’. (C'**)" (mod 2) (2) 


The formula below is transition matrix for the reverse of An, given thatn = 5k: 
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Msi? B B © B 


BT Mz! B 
M = BF AT No Noi p 
o Mz! B 
BT ok BT BT Mz? 
where 
00110 00110 
00011 00011 
Mz1ļ1 0 1 0 1|mod2) B=|00101 
11000 00000 
01100 00000 


|M, [mod2, representing the transition matrix starts with five numbers are shown below in (3). 
IM, |mod2 = {* ifn = 5k or n = 5k + 1,with k €E N) 
i 0, otherwise (3) 
2.2. Linear cellular automata transform 
This investigation involved doing the watermark embedding process on the vector map, which 
included the transformation frame of the coordinates of the vertices. This procedure was used on 
the coefficient of the transformation result frequency of the vector map’s data. Figure 1, shows that 
the coordinates are converted into an LCA transform, so that the vector map can be translated into 
the frequency domain of a signal. This is further supported by (4), which shows that the LCA transform, can 
convert the host map’s v,, coordinate. 


Watermark 
Data 


= LCAT Embedding iLCAT 
map 
R 5 


Embeded 
map 


Watermark LCAT — | 
Data 


Figure 1. Visualization of LCA transform 





< 





T(M) =Ni My. Vx1 (mod 2) (4) 


where, T(M)=host map's domain transformation value, M, = LCA's transition matrix, v,, = host map^ s 
digital media value, N = number of vertices altered by the frequency domain. 

Note: in (5) shows that the v,,; co-ordinate undertakes transformation with the LCA transform, allowing for 
the encrypted watermark part to be put into the equation N,N , vx, , TM), M,. 


Ver =Vx1' + aw (5) 


where, o=embedding parameter, W=watermark part. 

Variations of the vector map are directly proportional to the embedding parameter (a), meanwhile 
the watermark’s (W) resistances rises. This equation uses acceptable alterations to the vector map, a huge 
resistance value and 3-part a values. The inverse of the LCA transform is given in the following (6): 

In 6 below, represents the LCA transform’s inverse: v,, ” , iT(M) 
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iT(M) = XN- Mat. vy.” (mod 2) (6) 


where: iT(M) = host map's inverse domain transformation value , v,, "= digital media value of thehost 
map's transformation. 


2.3. Watermark embedding phase 

The suggested scheme is applied in the frequency domain and the digital watermark's embedded 
model is represented in Figure 2. The watermark insertion process consists of three steps, when a public key 
is used to encrypt (three parts, namely a vector map, the size of LCA transition matrix (Mn) and a watermark 
that scrambles the elements) and when a private key is used to decrypt. The proposed scheme relies on the 
linear cellular automata transform (LCAT) algorithm, and it operates in the frequency domain. Additionally, 
the proposed scheme does not fall under geometric attacks, namely, translation, scaling, and rotation. The 
algorithm of the scheme is as follows: 

Figure 2 displays the recommended scheme, used in the digital watermark’s embedded model and 
the frequency domain. There are three key steps to watermark insertion; when a public key is used to encrypt 
namely the size of Mn, watermark that mixes up elements and a vector map. This can happen when a private 
key is used to decrypt. This program operates in the frequency domain and depends on the LCAT (linear 
cellular automata transform) algorithm. The program is also is not susceptible to RST (Rotation, Scaling and 
Translation), which are considered to be geometric attacks. The algorithm goes like this: 

— Select two reference vertices, Vs, and vg, in the range (1 < Vf ,Vp2 < n ) as the vector map of M, in 
order to assure security. 

— When length (N) is transformed into a domain frequency (with no references), decide on the number of 
vertices in the map file (M). 

— Co-ordinates fare taken from the vertex and converted into a LCAT transform 

— Using equation 25, encrypt the factors of W*, which produces the data sequence 
W* ={w, Iw; € {0,1}, i=0, 1,...., 1—1} 

— Put W* into the final two consecutive digits (reduces the impact on the precision) and assume that a 
double floating point 16-digit coordinate value in a decimal fractional version. The embedded value is 
in the range of 0 to 99 and does not correlate with w;*. If we assume that the integer D is composed of 
the two digits, then: 


w= f w, is 0 then D < 50 and saved at the positions; } (7) 
7 w, = 1, otherwise 
— The initial shape file (with the inverse LCAT) of the frequency domain vector map is reinstated after 
the watermark is placed. 


2.4. Watermark extracting phase 
The watermark extraction point and the insertion process have similar steps, but in the opposing 
sequence. The three steps used in the exertion process are the results of the insertion process, namely the two 
reference vertices vf1 and vf2, fixed size LCA transition matrix (Mn), and the watermarked vector map, as 
shown in Figure 3. In the opposite sequence, the watermark insertion procedure and extraction point have 
comparable phases. The outcome of the insertion procedure, produces three steps for the exertion procedure. 
The results are Mn (fixed size LCA transition matrix), watermarked vector map Figure 3 and two reference 
vertices (vf1, vf2). 
— Select vg, and Vp (1 < Ve, ,Vep <n) , under private key k’s control, as the reference vertices for 
the vector map of M 
— Determine the number of vertices in the map file M and the length (N) that will later be transformed 
into a domain frequency without the references 
— Once the set of coordinates for each feature is obtained, it will be transformed into a LCAT transform 
— The embedded watermark location and watermark bits are extracted using (8): 
w= if D <50 then w,;,’ is 0 
7 | w;į* = 1, otherwise (8) 
— Extract the original embedded watermark sequence W, with private key k, by inversing 
the watermark pattern. 
— Get the watermark by rebuilding the watermark pattern. 
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— Select two reference vertices, Vs, and vg, in the range (1 < Vf ,Vf2 < n ) as the vector map of M, in 
order to assure security. 

— When length (N) is transformed into a domain frequency (with no references), decide on the number of 
vertices in the map file (M). 

— Co-ordinates of each feature are gained and then taken from the vertex and converted into 
a LCAT transform. 

— The (8) is used to extract the embedded watermarked parts and location. 

w= f D <50 then w,* is | 
wį* = 1, otherwise (8) 

— Reverse the watermark pattern, so that you can use the private key (k) to extract the original embedded 
watermark sequence (W). 

— Reconstruct the watermark pattern, in order to get the watermark. 


Original vector map Watermaked vector map 


(1) Extraction coordinates | 
of vertices 


(1) Extraction coordinates 
of vertices 


(Calculation of 
Matrix LCA 






(3) LCAT Transformation al 


Watermark 


(4) Eneryption of the 
Watermark 


(4) Extraction of Watermark 


(GQ) iLCAT Transformation 


(5) Generation of list bytes of 
Encryption Watermark 


(7) Generation Coordinates 
of Modified 


l 


Watermarked (6) Decryption of the 
vector map Watermark 





Figure 2. Watermark embedding Figure 3. Watermark extraction 


3. RESULTS AND DISCUSSION 
3.1. Experimental results 

For the evaluation of data type (.shp), there are two shapefile maps. Two vector maps make up 
the file type ESRI standard [32], which are a spot height and coastline map of Taylor Rookery [33]. 
Copyright marker was represented by a bitmap picture. The PC used for this study was Windows 10 
professional, QGIS version 3.0, 16 GB memory and 2.3 GHz. Programs used were python and MATLAB. 
The secret parts associated with each transform coordinate had a M, = 30, a in LSB and T = 1 for iterative 
embedding. The invisibility of our approach was displayed in the first test. Figure 4 shows vector maps 
watermarked, using the method suggested in section 3.5. This produced the watermarked types displayed in 
Figure 5. 

The recommended performance technique analysis for this study was measured by the NC 
calculation. It was also used to examine the similarities between the original watermark, before and after 
the extraction (values ranging 0-1). The success of the watermarking usage is indicated by how high the NC 
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value is. The larger the NC value is, the more similarity there is between images. In (9), w (initial value) and 
w;, are the extracted watermarks. Table 1 shows that the original and extracted watermarks are the same, as 
the evaluation result for NC is approximately 1, plus they have the same watermark content and length. 
This procedure was able to input copyrights as a watermark (without lessening its quality), as a watermark 
extracted again from the vector map file, does not need to go through dimensions or content changes. 


(a) (b) 


Figure 4. Test 2D vector maps: (a) spot height map of Taylor Rookery, (b) coastline map of Taylor Rookery 


(a) (b) 


Figure 5. The watermarked 2D vector maps of Figure 4 map 


M * 
di=0 Wi X Wi 


foxtocwye x, fowo 0) 


Table 1. Result of similarity test between original watermarks and extracted watermark 
Map Map Type Features/vertices Original Watermark Extracted Watermark NC 


À 
À 


NC = 


ae 





a Point 355/355 0.998909 


b Polyline 18/4279 0.998308 





3.2. Invisibility evaluation 
Table 2 shows the two parameters that were used previously in the invisibility measurement 
(as reference analysis for calculating RMSE), which subsequently decides on the alteration between 
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the interpolated watermark’s results and the beginning of the map file. The equation for the making of 
the RMSE, is shown below: N,v ,,*”" 


1 * TT 
RMSE =(— Lito vivit”) (10) 


where, N =the number of vertex map vectors, V,, =equivalent x coordinates in the initial vector map. 


* TI 


V y1 «= equivalent x coordinates in the recovered vector map. 


Table 2. RMSEs between the recovered maps and the original maps (x10? (T = 1) 


Map DCT DWT FFT Proposed (LCAT) 
a 5.0917 10-4 1.3014 10-4 2.9023 10-5 1.1317x 10-9 
b 1.8514 10-3 3.5872 10-4 2.3342 10-5 3.1562 10-9 


3.3. Fidelity evaluation 

Digital watermarking is untraceable by human senses, making it a very reliable system, plus it 
doesn’t significantly worsen the media file interpolation. Moreover, with regards to RMSE, the furthest 
changes can be gauged. A position change is indicated by the furthest distance, which is caused by 
watermark interpolation into vector files. This is gained via an assessment comparing the original vector map 
file, vector map file (containing the watermark) and coordinated vertex. To convert the furthest distance into 
meters, QGIS is implemented. The lengthiest noticeable shift in data analysis is 60 cm, as seen in Table 3, 
therefore the furthest distance value still maintains a significance level of precision in the vector map, which 
is displayed in Figure 6. 


Table 3. Result of fidelity test between original maps with watermarked mas 


Map T (m) DCT DWT FFT Farthest Distance (meter) LCAT 
a 0.5 0.45 0.4 0.3 0.050 
b 0.5 0.5 0.43 0.38 0.068 
TE ee — 
= r : ¥ Py i : i E i “ i 
(a) (b) 


Figure 6. Original and watermarked map overlaid with each other 


3.4. Robustness evaluation 

Examining the strength of the watermark is the next stage. Advanced methods are used to challenge 
the three geometric attacks, classed as RST (rotation, scaling and translation). The RST also includes a signal 
operation attack of three components; vertex deletion, vertex modification and vertex insertion. Attacks in 
the form of a vertex point embedded in watermarks and vector maps with spatial features, were used to test 
on the data. In order to conduct attacks on the test vector maps, QGIS software is used and 
the results of the extraction to determine the impact of the attacks Table 4, are found out by doing NC 
calculations. Table 4 can be distinguished into three parts: 1) Rotation attack: rotating the test data 
coordinates by 30° to 180°, 2) scaling attack: enlarging the test map from 0.25 to 4.0, and 3) translation 
attack: Alter the positioning of the vertices by moving the coordinates. 
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Table 4. Experiment results for robustness 


Attack Type Detection Rate (NC) 
Rotation (p=30°) 0.939542 
Rotation (p=60°) 0.921861 
Rotation (p=90°) 0.906505 

Rotation (p=180°) 0.898085 
Scaling (s = 0.25) 0.913221 
Scaling (c = 0.5) 0.890227 
Scaling (c = 2.0) 0.884843 
Scaling (c = 4.0) 0.880051 
Translation (—1.2 m, -2.3 m) 0.901241 
Translation (4.2 m , 5.6 m) 0.896940 
Translation (2.6 m, — 7.9 m) 0.765425 
Translation (— 0.6 m, 0.7 m) 0.558168 


3.4.1. Rotation attack 

Turning a vector map, in order to reduce the WDR (Watermark Deduction Rate) is the definition of 
rotation. Figure 7 shows that there was a rotation from p=30° (low point) to p=180° (high point). Only very 
small effects on the quality of the embedded watermark were shown by the testing. We now have evidence 
that the recommended WV map is sufficient enough to resist such an attack. 


WATERMARK DETECTION RATE ON ROTATION ATTACK 


Da - 0.939542 
0.94 
0.935 - 
0.93 


0.925 - 





0.92 


0.915 - 







0.91 - “906505 


DETECTION RATE (NC) 


0,905 


, 0.858085 
0.9 - 








0,895 


* 


n (p=180°) 


0 50 100 7 209 Rotation (p=90°) Rotatio 


Rotation 


Figure 7. Quality of watermark vector map after rotation attacks 


3.4.2. Translation attack 

Scaling (for this study), is defined as resizing a map along its axes, so that the WDR can be reduced. 
The scaling range was ¢ = 0.25 (low) to c = 0.25 (high). Akin to the rotation attacks, scaling attacks also have 
only a small impact on the WDR and Figure 8 shows that the watermark embedded, just about managing to 
survive. We can conclude that under a scaling attack, this method of LCAT watermarking is only 
modestly robust. 


3.4.3. Translation attack 

Relocating the whole map shifting it into a specific direction, is the main trait of the translation 
attacks. The map used in this research, translated from -1.2, -2.3 m to -0.6 m, 0.7 m in Figure 9. The studies 
only showed a minor effect on the quality of embedded watermark, apart from the exceptions at -7.9, -0.6 m, 
0.7 m, and 2.6. The minimal impact could be due to the translation attacks, so we can conclude that 
the LCAT watermarking technique is moderately resistant to translation attacks. Watermark extraction was 
not always attained, by the attack techniques (with NC=1). Out of the entire evaluation scenario, the failure 
rate of the watermark detection was 89%. From this, we know that these extraction methods can’t be used in 
attacks that change the watermark bit (Table 4). Friction against the vertex is made by TAs, changing 
the WBV (Watermark Bit Value) and stopping extraction. In the meantime, rotation attack results show that 
the watermark extraction will be unsuccessful is the bit value changes by up to 1. If the LCAT algorithm 
can’t maintain its fixed genuine watermark level during an attack, then the WBV changes and extraction 
becomes unattainable. 
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WATERMARK DETECTION RATE ON SCALING ATTACK 














0.915 „0-913221 

0.91 - 
g 0.905 
cs 
w 0.9 - 
z 0.895 - 
O 
5 0.89 - 
F 0.884843 
A 0.885 - 

0.880051 
0.88 
0.875, +a 
0 0.5 1 1.5 2 2.5 3 3.5 4 
Scaling 
Figure 8. Quality of watermark image after scaling attacks 
WATERMARK DETECTION RATE ON TRANSLATIN ATTACK 
i E er pam ee TRANSLATION 
0.9 
0.76542 


.558168 


DETECTION RATE (NC) 
oo 
[S] 


a] 
à 





(—1.2 m, -2.3m) (4.2,5.6m) (2.6,—7.9) (-0.6 m, -0.7m) 


Figure 9. Quality of watermark vector map after translation attacks 


3.5. Comparative study of watermarking techniques 

There are various watermarking techniques describes in this study. When they are tested on DWT, 
FFT and DCT algorithms; fidelity, geometrical attacks and better invisibility is shown by the linear cellular 
automata [31]. Following this, 1D image files were watermarked by CAT; despite the fact that 2-D CA 
needed to be converted between two dimensions, specifically along with repetition. The same keys, vertices 
and maps shown in Table 5 were used, to compare them to our program [31]. Table 6 summarizes 
the resistance of the four techniques against the RST attacks (in relation to watermark immunity to 
elimination and inadvertent degradation). LCAT comes out on top, as a detection rate of 0.76; indicates that a 
complete watermark recovery is still feasible. 


Table 5. Parameters analysis of both DCT, DWT, FFT and LCAT 


Parameters analysis 


Evaluation Metrics LCAT DWT FFT DCT 
Invisibility high medium high medium 
Fidelity high Low medium Low 
Geometrical Attacks 89% 717% 83% 174% 
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Table 6. Watermark detection rate comparison for DCT, DWT, FFT, LCAT under various attacks 


Attack Type DCT DWT FFT LCAT 
Rotation (p=30°) 0.87 0.85 0.88 0.94 
Rotation (p=60°) 0.84 0.86 0.90 0.92 
Rotation (p=90°) 0.80 0.80 0.86 0.91 

Rotation (p=180°) 0.62 0.80 0.79 0.90 
Scaling (s = 0.25) 0.47 0.62 0.72 0.91 
Scaling (c = 0.5) 0.65 0.70 0.78 0.89 
Scaling (s = 2.0) 0.87 0.90 0.93 0.88 
Scaling (c = 4.0) 0.78 0.73 0.84 0.88 
Translation (—1.2 m, -2.3 m) 0.77 0.83 0.88 0.90 
Translation (4.2 , 5.6) 0.82 0.83 0.87 0.90 
Translation (2.6, — 7.9) 0.71 0.73 0.80 0.85 
Translation (-0.6 m, -0.7 m) 0.68 0.60 0.67 0.76 


Manipulation of results on the watermark insertion map, is the cause of watermark extraction 
failure. The distortion value is affected by the use of the limitation value in the watermark insertion. 
The watermark value is found by using a bit matrix size of M, < 30, in order to allow for extraction. 
Limitation value is affected by modification amplitudes of M, = 35 bits. The techniques of some of the tests, 
managed to keep the same watermark, extracted and original watermarks were the same and NC = 1, in spite 
of the low robustness. For every LCAT value, modifications on the value of the sequence complex on 
the vector mapping, an LCAT calculation was spread for each one. This was done to keep the LCAT value 
within the limits of the watermark extraction. Alterations to the coordinate value, directly affects the inserted 
watermark bit value; when insertion happens on the spaced-out domain. This will produce contrasting results. 
The transform domain calculation spanning the LCAT is better at conserving the watermark, in comparison 
to the spatial technique. The robustness of the method depends on, FDA (frequency domain algorithm), 
extraction limit, related programming methods, asymmetric algorithm key quality and the data storage length. 


4. CONCLUSION 

Finding a good frequency watermarking scheme for 2D vector map drawings was the principal aim 
of this paper. The frequent problems with watermarking schemes of vector maps are invisibility, high fidelity 
and weak robustness. Furthermore, the scheme needs the original cover when extracting the watermark, 
because it is ‘Non-blind’. A watermarking technique (that was LCAT domain transformed) was embedded 
into the vector map. Invisibility showed that likeness in fidelity stages, in the watermarked map are upheld. 
The RMSE (distortion scale) staying circa zero and the distance from the original vector map being 10% or 
less, also displayed this. The impact of the frequency domain scatters throughout various frequencies, 
restraining impact and upping this method’s dependability. LCAT FDA warrants the vector map’s exactness 
and robustness is shown in 89% of examined cases. 
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